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Abstract 

We solve a 4-(bond)-vertex model on an ensemble of 3-regular ($ 3 ) planar random graphs, which 
has the effect of coupling the vertex model to 2D quantum gravity. The method of solution, by 
mapping onto an Ising model in field, is inspired by the solution by Wu et.al. of the regular lattice 
equivalent - a symmetric 8-vertex model on the honeycomb lattice, and also applies to higher valency 
bond vertex models on random graphs when the vertex weights depend only on bond numbers and 
not cyclic ordering (the so-called symmetric vertex models). 

The relations between the vertex weights and Ising model parameters in the 4- vertex model on 
$ 3 graphs turn out to be identical to those of the honeycomb lattice model, as is the form of the 
equation of the Ising critical locus for the vertex weights. A symmetry of the partition function 
under transformations of the vertex weights, which is fundamental to the solution in both cases, can 
be understood in the random graph case as a change of integration variable in the matrix integral 
used to define the model. 

Finally, we note that vertex models, such as that discussed in this paper, may have a role to play 
in the discretisation of Lorentzian metric quantum gravity in two dimensions. 



1 Introduction 



Random matrix models have proved remarkably successful in investigating the critical behaviour of various 
sorts of spin models (Ising ||, g], Potts ||], O(N) [§)■■■) living on ensembles of planar random graphs. 
This, in effect, couples the models to 2D quantum gravity which originally motivated their study in the 
context of string theory and random surfaces H . From the purely statistical mechanical point of view 
one is looking at the critical behaviour of an (annealed) ensemble of random graphs decorated by the 
spins. The connection with gravity comes from the fact that the graphs can adapt their connectivity 
in response to the spin configuration, which in turn is influenced by the connectivity of the graph on 
which it lives. This captures the essential feature of the back reaction of the matter on the geometry that 
characterises gravity in a discrete form. 

The graphs of interest appear as planar Feynman diagrams in the perturbation expansion in the vertex 
coupling of a Hermitian matrix model free energy of the general form |j| 

F= ^ iog /n^ ex p(- 5 ) w 

i 

in the limit M — > oo, where the 3?j are sufficient M x M Hermitian matrix variables to describe the 
matter decoration and the log ensures that only connected graphs contribute. For the Ising model on 
trivalent planar graphs, for instance, we need two matrices X, Y and the action takes the form 



S = Tr^(X 2 + Y 2 )-gXY-^ [e h X 3 + e^'F 3 ] , J 



(2) 



where g = exp(— 2/3) and h is the external field. We can see that the X 3 terms can be thought of as 
spin- up or + vertices and the Y 3 terms spin-down or — , while the inverse of the quadratic terms gives 

the appropriate edge weights for + +, H — and edges. 

To date there has been less consideration of arrow vertex models on ensembles of planar random 
graphs and none at all of bond vertex models. In the arrow case one decorates the edges of the graph 
with arrows and assigns weights to the vertices depending on the local arrow configurations around each 
vertex, in the bond case one simply has empty or occupied edges (bonds) and vertex weights depending on 
the local bond configuration. Some arrow vertex models have been discussed using elementary means in 
(?|, mostly on 4-valcnt random graphs, revealing XY and Ising critical behaviour for suitable coupling- 
constant loci. More recently some elegant and sophisticated work using character expansions || has 
provided a full solution of an arrow vertex model on 4-valent planar graphs confirming the XY critical 
point found in 

In this paper we give the first solutions of bond vertex models on planar random graphs. Our main 
example is a 4-vertex model on $ 3 random graphs, but the general method works for any bond vertex 
models whose weights depend only on the number of bonds at a vertex rather than their particular 
cyclic ordering. On regular lattices such models are termed symmetric vertex models and we retain the 
nomenclature here. Since our 4-vertex model on <3> 3 random graphs is the random graph equivalent of the 
symmetric 8-vertex model on the honeycomb lattice discussed by Wu et.al. [|To|, |ll|, [l^] we take our cue 
from where the honeycomb lattice model was first solved by making use of the so-called generalised 
weak graph transformation jl3| between vertex weights. This brought the partition function to a form 
that was recognisable as that of an Ising model in field an the critical behaviour of the model was then 
deduced from the corresponding Ising results. 

In what follows we first define our 4-vertex model, then solve it by mapping onto the Ising model in 
field on <I> 3 planar random graphs, which itself was solved in We move on to discuss higher valency 
symmetric models and show that the same general method of solution still works. We then compare 
these results with those of |T(], O |l2) for the 8-vertex model on the honeycomb lattice and note that the 
generalized weak graph symmetry of the honeycomb lattice model finds a counterpart in an orthogonal 
variable transformation in the matrix integrals used to define the random graph models. Finally, we close 
with some speculations on the relations of vertex models to Lorentzian metric quantum gravity, where 
they can be thought of as providing a discretisation of causal structure. 
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2 The 4- Vert ex Model on 3-valent random graphs 



The canonical (fixed number of vertex) partition function for the 4-vertex model on <I> 3 planar random 
graphs is given by 

Z 4v = Z 4v (a, b,c,d)=J2J2 a Nl t> N *c Ns d Ni (3) 

$N G 

where there are two summations: over <E> 3 random graphs with N vertices and ^ G over all possible 

bond configurations on each <E> 3 graph built using the vertices of Fig.l. Each graph has N\ vertices of 
type (a), N 2 of type (b), N 3 of type (c), N A of type (d) and N = Nt + N 2 + N 3 + 7V 4 . 

We can generate the graphs we require, suitably decorated with the allowed bond configurations, by 
inserting the action 



S' : : Tl \ ^(X 2 +Y 2 ) - A 



^-X 3 + ^Y 3 + cXY 2 + bX 2 Y 



(4) 



into equ. (|I]) , where X, Y are M x M Hermitian matrices and we have introduced an overall vertex coupling 
A for convenience.. The limit M — > oo picks out the planar diagrams. A direct attempt at solving a 
matrix model with the action equ. (0) might be possible using the character expansion techniques of H , 
but it turns out to be easier to follow in the footsteps of |n], |ll[ and establish a correspondence with an 
Ising model in order to determine the critical behaviour of the vertex model. 
We start by by carrying out the orthogonal transformations 



X -» (X + Y)/V2 

Y -> {X-Y)/V2 (5) 



on the matrices in the Ising action of equ.(||), followed by the rescalings X — > X/(l — g) 1 ^ 2 , Y — > 
Y/(l + g) 1 / 2 , A — > y/2X(l — g) 3 ^ 2 - This change of variable in the matrix integral has a trivial Jacobian 
and gives the new action 



[1, , 9x Acosh(/i) r , „ ,-. \smh(h)(q*) 3 / 2 
S = Tr \-(X 2 + Y 2 ) -^-[X 3 + 2,g*XY 2 }- 



Y 3 + — X 2 Y 



(6) 



where g* = (1 — + .<?), which is clearly of the same general form as equ.(Q). The bond vertex 

weights from equ.(|6[) are also given in Fig.l. It should be remarked that the transformation of equ.(||) 
is the duality transformation for the Ising model on random <I> 3 graphs JlJ], and that the bond graphs 
generated by the transformed action in equ.(^) are precisely the high-temperature expansion graphs of 
the original Ising model on $ 3 graphs. 

We can see that the vertex weights coming from the Ising model in equ.(|6|), which we denote by tilde'd 
quantities, 

a = cosh(/i) 

b = siiMh)(g*) 1/2 

c = cosh(h)g* 

d = sinh(»( 5 *) 3/2 (7) 

are fixed in the ratio ad — be, which will not be the case for the generic weights of equ.(|4|). However, if 
we consider the effect of the further orthogonal rotation Q 



(8) 



X \ f cos(<9) sin(0) \ / X 
Y ) ^ \ -s\n{6) cos(6>) )\Y 

^^This preserves the propagator X 2 + Y 2 in the action of equ.(^) and again has a trivial Jacobian. so it does not effect 
the evaluation of the matrix integral. 
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on the vertex weights of equ.M) 



we find 
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Ky 2)3/2 



a + 3yb + 3y 2 c + y 3 d] 
-ya+(l~2y 2 )b~(y 3 -2y)c- 
(y 3 ~2y)b+(l-2y 2 )c~i 
f 3y 2 b - 3yc + d] 



\v 2 a 



-y 2 d] 
yd] 



(9) 



-j/V 



where we have extracted cos(#), denoted tan(#) as y and deliberately used tilde'd quantities again for 
the transformed weights. We can now demand that these transformed weights satisfy ad = be, just as in 
equ.(|6|) which has the effect of making the partition function of the vertex model equivalent to that of 
the Ising model in field. Setting ad = be gives the following equation for the rotation parameter y from 
equating the right hand sides of equ. (9) 



By 2 + 2{C - A)y - B = 



(10) 



where A = bd — c 2 , B = ad — be, C 



ac 



b 2 . Since — oo < y < oo is at our disposal we can solve this 



for general a, b, c, d to obtain transformed weights satisfying ad = be. 

To deduce the critical behaviour of the vertex model by using this correspondence with the Ising 
model in field we need the relation between the vertex weights and the Ising parameters f3 and h. The 
Ising and tilde'd vertex model parameters were related by 



tanh(/3) 
tanh(/i) 



(11) 



as can be seen directly from equ.(|] 
original weights and y using equ.( 9 
Ising parameters (3, h 



We can then substitute for the tilde'd weights in terms of the 
to give the relation between the original weights a,b,c,d and the 



exp(2/3) 
tanh(/i) 



By + C-A 



A + C 



W 



By + 2C 
By-2A 



1/2 



(12) 
a + c and 



where y is again the one of the solutions of equ. (10), A, B,C are as above, T — (b + d)y 
W =(b + d)-{a + c)y. 

On <I> 3 random graphs the Ising model displays a third order magnetisation transition in zero field. 
The second of equs.(|l2|) shows that zero Ising field implies that W — 0, which in terms of the vertex 
weights (when y is substituted for) gives the locus 



a(b 3 + d 3 ) - d(a 3 + c 3 ) + 3(ab + be + cd){c 2 - bd - b 2 + ac) = 0. 



(13) 



From the vertex model perspective it is most natural to disallow loops of length one and two in the 
random graphs in order to make sure that there are no double edges or tadpoles, in which case the Ising 
critical temperature is given by f3 c — \ log (^) = 0.7733185 (2). The first of cqus.(^2|) then gives 



103 _ y/(ac -b 2 -bd + c 2 ) 2 + (ad - be) 2 
~28" ~ ac-b 2 + bd- c 2 



(14) 



as the position of the third order critical point on the locus equ.([13|). We still have, of course, the echo 
of the field driven Ising transition appearing for (3 > (3 C along the locus equ.(fL3|). 
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3 Higher valency vertex models 



In this section we show that on 4-valent ($ 4 ) random graphs the Ising equivalency can still be easily 
established by rotating the variables in a matrix model action. One extra ingredient in the 4-valent case 
compared with the 3-valent case is that the cyclic order of the bonds allows one to distinguish between two 
different sorts of vertices with two occupied and two empty edges, one with occupied bonds at right angles 
and one with occupied bonds "straight through" . The Ising equivalence only holds for the symmetric 
case in which equal weights are assigned to these bond configurations. 
On <I> 4 random graphs the Ising model action is 



A 



S = Tr I -{X 2 + Y 2 ) - gXY + - [e h X 4 + e^F 4 ] , 



(15) 



with only a higher order potential distinguishing the action from that for $ 3 graphs in equ. (j^) . The change 
of variables X -> (X + Y)/V%, Y -> (X - Y)/V2 and rescalings X -> X/(l - g) 1 ' 2 , Y -> Y/(l + g) 1 / 2 , 
A ^ 2A(1 — g) 2 gives the action 



S = Tr\^-{X 2 + Y 2 ) 



XC0 f {k) [X 4 + [g*) 2 Y 4 ] + ACQSh o Wg * [2X 2 Y 2 + XYXY] 



+ Xsimh(h)g* 



g*XY 



X 3 Y 



(16) 



where, as on $ 3 graphs, g* = (1 — g)/(l + g). Note the presence of the [2X 2 Y 2 + XYXY] term fixing 
the weights of right angled and straight through two-bond vertices to be the same. It should be remarked 
once again that the change of variables in equ.(|^) is a duality transformation for the model, and that 
equ. (16) generates the high temperature expansion graphs of the original Ising model in equ. (|l5|) . 

The effect of the orthogonal rotation of equ.(^) on the vertices arranged in some suitable lexicographic 
order (e.g. X 4 , X 3 Y, X 2 Y 2 , XY 3 , Y 4 as labelled on Fig.2) gives 

[a + Ayb + 6y 2 c + Ay 3 d + y 4 e] 
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[-ya + (1 - 3y 2 )b - (3y 3 - 3y)c + (3y 2 - y 4 )d + y 3 e] 
[y 2 a + (2y 3 - 2y)b + (y 4 - Ay 2 + l)c + (2y - 2y 3 )d + y 2 e] 



-y 3 a + (3y 2 - y 4 )b - (3y - 3y 3 )c + (1 - 3y 2 )d + ye] 



[y 4 a - 4y 3 b + 6y 2 c - 4yd + e] 



(17) 



(18) 



The strategy for solving the $ 4 model is identical to the $ 3 case. One uses the above transformations 
of the vertex weights to bring the model to the Ising action defined above in cqu.(|l6|). The tilde'd weights 
in this (a for X 4 etc.) can be characterised in various ways, for example: ad — dc, c 2 — de. Having 
achieved this the critical behaviour can again be read off from that of the Ising model in field, this time 
on $ 4 graphs. Indeed, it is clear that generic symmetric vertex models on random graphs can be solved 
by using this method. In the q-valent case one carries out the same sequence of two orthogonal rotations 
on the action for the Ising model on <I> 9 planar random graphs 



S = Tr l^(X 2 + Y 2 ) - gXY + ^ [e h X q + e - h Y q ] J 



(19) 



to obtain the action for the symmetric vertex model. One can thus deduce that the critical behaviour of 
the symmetric vertex model on q-valent random graphs will be Ising-like. 
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4 Comparison with the Symmetric 8- Vertex Model on the hon- 
eycomb lattice 

The closest regular lattice equivalent of our <I> 3 planar random graphs is the honeycomb lattice where each 
vertex still has coordination number three but, in contrast to the random graphs, all loops are of length 
six. The eight possible bond vertex configurations all look similar to those for the random lattice in Fig.l 
but we now have orientational order as well as cyclic order around each vertex so b and c represent three 
possible orientations each. However if we impose equality between the different orientations, we arrive at 
the symmetric vertex model of jlO[ [TlJ [l2| which, like the $ 3 random graph vertex model, has 4 distinct 
vertex weights. The partition function is given in a similar fashion to equ.(^) by 

Z Sv = Z Sv (a, b, c,d) = Y^ a Nl b N2 c N3 d N * (20) 

G 

where there is now no sum over random graphs since all the bond configurations are living on the 
honeycomb lattice. 

Several symmetry properties of the partition function are immediately apparent. From the nega- 
tion of vertices occurring in pairs we have Z% v {a,b,c,d) = Zs v (—a,—b,—c,—d) = Z$ v (—a,b,—c,d) — 
Z$ v (a,—b,c,—d). Similarly, exchanging dark and light bonds gives Z% v {a,b,c,d) = Z% v {d,c,b,a). Per- 
haps rather less obvious is the generalised weak graph symmetry JlJ] Zg v (a,b,c,d) = Zg, v (d,b,c,d), 
where the transformations between the original and tilde'd variables are precisely those generated by the 
orthogonal rotations in the random graph model in equ.(||). 

In the generalised weak graph symmetry on the honeycomb lattice was thought of as being 
generated by rotations with 

v(y) = „ , * ( ! \ ) (2i) 



(i + y 2)i/ 2 y y -i 

on the vertices in "bond space" , where occupied and empty bonds on a given edge were to be thought of 
as components of a vector. An alternative choice of transformation 

U(V) = T, , !.2M/ 2 ( L ? ) ( 22 ) 



(l +1/ 2)l/2 ^ -y 1 

was related to V(y) by the negations b — > —b, d — > — d and was hence not independent. Remarkably, the 
second transformation is just that generated by equ.(||) in the random graph model. 

The generalised weak graph symmetry played an important role in the original solution of the hon- 
eycomb lattice model pd[ because a suitable choice of y allowed generic vertex weights a, 6, c, d to be 
transformed to a, 6, c, d satisfying ad = be just as for the random lattice. For vertex weights satisfying 
the latter condition 

Z 8v (a~b,c,d) = (~acosh(h)/2) N ( cosh(P))- 3N / 2 Z Ismg ((3,h) (23) 

where Zj sing was the standard honeycomb lattice Ising partition function. The equivalence followed from 
observing that this choice of weights meant that the vertices generated precisely the right weights for the 
diagrams of the high temperature expansion of the Ising model on the honeycomb lattice. 

The critical behaviour of the symmetric 8- vertex model on the honeycomb lattice was thus deduced by 
following exactly the same path we have taken for the 4-vertex model on $ 3 random graphs - exploiting 
a symmetry of the partition function to transform it to that of an Ising model in field. The parallels run 
even closer since equs.([| [To|, |l2] , |l3| ) still hold identically for the honeycomb lattice vertex model and 
the only difference in equ.(|l4 ) is in the numerical value of the Ising critical temperature on the left hand 
side. This might grandiloquently be phrased as a non-renormalization theorem for the vertex weights 
since the ratios of weights in equ.(|ll|) are unaffected by the coupling to gravity that is represented by the 
sum over random graphs. 

The close correspondence between the ( 1 )3 and honeycomb results becomes a little less surprising when 
one considers a second determination of the critical behaviour of the honeycomb lattice model by Wu in 
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where it was shown that a direct mapping between the Ising and vertex models was possible when a 
decoration- iteration transformation was employed. In Fig. 3 we show a single vertex, which hosts a spin 
subject to an external field H . In addition the edges have spins which are subject to a field 2H' and the 
two sorts of spins interact via a spin-spin coupling R. The presence of a dark bond can be denoted by an 



edge spin a — 1 and its absence by a 
in terms of appropriate H' , H and R 



-1. The vertex model weights can then be represented correctly 



cosh(2i?) 
exp(4ff') 
cosh(2if) 




AC \4AC 



B 
46c" 



(24) 



when the central vertex spins are traced over. One can go in the other direction and decimate the edge 
Ising spins by replacing the two R interactions and edge field 2H' with a single interaction (3 and fields 
at the end of each edge (i.e. on the vertices). This gives one a standard honeycomb lattice Ising model 
in field. The relation between the vertex model weights and the Ising parameters f3, h determined in this 
manner is identical to that cqu.( 12). 

This decoration-iteration approach strongly supports the observation here that the equivalence be- 
tween the symmetric 8-vertex model on the honeycomb lattice and the Ising model in field continues to 
hold for the 4-vertex model on <I> 3 random graphs. The decoration iteration transformation is a local 
transformation which relies only on the valency of each vertex and in both cases all vertices have valency 
three. It also supports the result that the relation between the Ising and vertex parameters is identical 
on $ 3 random graphs and the honeycomb lattice, since the decoration-iteration transformation doesn't 
"see" the randomness, so as far as it is concerned the random and honeycomb lattices are identical 0. 

The critical behaviour of symmetric vertex models on higher valency regular lattices may be deduced 
by using the appropriate generalized weak graph transformations to map the models onto Ising models 
in field in all cases. The generalised weak graph transformation written down by Wu for general valency 

q 



(l + y 2 )i/'- 



E 

k=0 



q-i 
j -k 



-l) k y i+j - 



21; 



(25) 



where the vertices are labelled in some suitable order as in the $ 3 , $ 4 cases, is precisely that obtained by 
picking out the appropriate powers in an expansion of terms such as (cos(6)X + s'm(8)Y) q ~ k (— sin(6)X + 
cos(6»)F) fe . 



5 Speculations on Connections with Lorentzian Gravity 

$ 3 graphs are duals to triangulations which are the natural discretisation of Euclidean signature (+ +) 
2D spacetimes. An analytical continuation to physical Lorentzian (— +) signatures is by no means on 
such firm ground as the analogous procedure in standard quantum field theory on flat backgrounds. It is 
therefore of some interest to attempt to formulate discretised theories which might serve as toy models 
for Lorentzian gravity. 

An essential extra ingredient compared with the Euclidean case for any Lorentzian theory is a causal 
structure. One then faces the option of whether causal structures should be summed over or imposed. A 
relatively rigid choice was explored in p^ j where triangulations with spacelike slices joined by fluctuating 

2 One possible fly in the ointment is that loops of length one and two, i.e. double edges ("bubbles") and self-loops 
("tadpoles"), are in principle present in the random graphs. The star-triangle relation would degenerate at vertices adjacent 
to such edges. The simplest solution, which we have taken, is to restrict the ensemble of random graphs to exclude such 
configurations, which does not affect the critical behaviour of the Ising model though non-universal quantities such as the 
critical temperature may be changed. 
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timelike edges were used to triangulate Lorentzian space. It was found that only when branchings were 
allowed was the critical behaviour of Euclidean gravity recovered upon analytical continuation back. 

Another possibility is to postulate that global, or at least long-range, causal structure might emerge 
dynamically in much the same way that the fractal structure of the 2D manifolds does in the Euclidean 
case. Inspired by ||l6| , |l7f we introduce local causal structure by using spacelike edges to triangulate our 
Lorentzian spacetime. The normals to all the triangles are then time-like, which can be either past or 
future directed. It is then natural to return to the dual spin network picture, which gives an embedded 
directed <i> 3 graph if we denote the causal relations between points in the interior of adjacent triangles by 
arrows pointing from the past to the future of each vertex. This should be contrasted with Euclidean 2D 
triangulations where the dual picture is of undirected <i> 3 graphs. 

It is thus natural to use arrow vertex models in this context, rather than the bond vertex models 
considered so far. However the translation between arrows and bonds turns out to be trivial, as we see 
below. The various sorts of 3-valent (<I> 3 ) arrow vertices are most conveniently labelled by their (indegree, 
outdegree). The different possibilities are (3, 0), (0, 3), (2, 1) and (1, 2) as shown in Fig. 4. A natural choice 
is to take conjugate weights for the (3, 0) and (0, 3) vertices and similarly for the (2, 1) and (1, 2) vertices 
as this preserves the symmetry under arrow reversal. The partition function for this 4-arrow- vertex model 
is thus of the form 

z =Y,J2 aNl ^ N2bN3lNA ( 26 ) 

where the first sum is again over different planar <I> 3 graphs, the second over arrow assignments and 
we have N\ (3,0) vertices, N2 (0,3) vertices etc. Note that only the (2,1) and (1,2) vertices can be 
considered as regular spacetime points since the (3, 0) and (0, 3) vertices act as microscopic black and 
white holes respectively. 

We can obtain this partition function from the N vertex term in the expansion of the free energy of 
the complex matrix model with the action 

S = -Tr - - [($ f ) 3 + $ 3 ] - i- [($ 3 - ($t) 3 ] - 7 [(^) 2 $ + $ 2 ($t)] - iS [($t) 2 $ - $ 2 ($ t )] ) (27) 

2 13 3 J 

where is now an M x M complex matrix and a = a + i/3, b = 7 + iS. We can transform simply to 
a bond vertex formulation of the model by splitting $ into Hermitian components X + iY . This allows 
us to interpret Y edges as those containing bonds and X edges as empty (or vice-versa). The resulting 
action in terms of X, Y 

S = Tr {-{X 2 + Y 2 ) - (Q + 37) A 3 - (/? + 3(5) y 3 - (7 - a)XY 2 ~ (5 - (3)X 2 y) . (28) 
12 3 3 J 

is now clearly that of our 4-bond- vertex model on $ 3 random graphs. This displays Ising criticality, as 
we have seen, so a suitable tuning of the vertex couplings in equ. (p8f) (and hence in the original arrow 
vertex model of equ. (p7|) ) can reach this point ^|. Since this is a continuous transition one might hope 
that the diverging correlations of the Ising model at the critical point might be translated back to the 
appearance of long range casual structure ex nihilo in the arrow vertex model formulation. However, one 
can see heuristically that closed loops of arrows are present at all scales at the transition point. Cluster 
algorithms exist for vertex models Jl8| which act by identifying such closed loops and flipping them. 
These are effective at combatting critical slowing down so they are identifying and flipping loops at all 
scales near a continuous transition point. The loops are innocuous from the statistical mechanical point 
of view, but they represent closed timelike loops when the vertex arrows are interpreted as giving the 
causal structure. 

A profusion of closed timelike loops is obviously an undesirable property and we thus conclude that 
the particular vertex model discussed here is probably not a suitable candidate for modeling discretised 
Lorentzian gravity as it stands. It would be interesting to determine what modifications might make it 
so. 



interestingly, the presence of the black hole and white hole vertices appears to be necessary to do this, there is insufficient 
freedom with only the regular vertices. 
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6 Discussion 



We have solved a bond 4-vertex model on an ensemble of planar $ 3 graphs by taking our lead from the 
honeycomb lattice solution of the symmetric 8-vertex model and exploiting its equivalence to an Ising 
model in held. An important ingredient of the solution was an orthogonal rotation in the matrices use to 
dchnc the random graph model in order bring a generic weight configuration on to the Ising locus. This 
turned out to be functionally identical to the generalised weak graph transformation used in the original 
honeycomb lattice solution |hJ . The method of solution and the equivalence of the matrix rotation to 
the generalised weak graph transformation on an equivalent regular lattice were also shown to work for 
higher valency symmetric vertex models on random graphs. 

We have seen that the relations between the vertex model parameters and the Ising parameters were 
identical for the honeycomb and random graph models. This is not so surprising at it might first seem 
when viewed in the light of the decoration-iteration solution on the honeycomb lattice, since this depends 
only on the valency of the vertices, which is identical in the random case (modulo caveats about tadpoles 
and bubbles discussed previously). This result could be couched as a non-renormalization theorem for 
the appropriate ratios of vertex weights, since putting the models on ensembles of random graphs is 
equivalent to coupling them to 2D quantum gravity. 

The methods discussed in this paper are restricted to symmetric vertex models since the orthogonal 
rotation of equ. (||) necessarily gives only symmetric weights when applied to Ising actions such as those in 
equs.((^, |l5| , [l9| ). From the matrix model point of view the vertex model solution shows that an apparently 
hopeless potential containing X 3 ,Y 3 , X 2 Y and XY 2 terms can still give rise to a soluble model, once 
the matrices are appropriately transformed. Finally, we have speculated on the relation between vertex 
models and Lorentzian signature gravity, pointing out a potential problem with closed timelike loops 
when employing the class of models considered here. 
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a 



Figure 2: The 5 different vertex weights in the symmetric 16 vertex model on the square lattice. The 
"straight-through" two bond configuration which is not shown also receives weight c. 
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Figure 3: The decorating spins on the edges are subject to an external field 2H' , the central spin to a 
field H. The edge-vertex spin interactions have weight R. 
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(a) 0» 




(c) (d) 



Figure 4: The possible arrow vertices which appear in the model. Their weights are equivalent to linear 
combinations of the bond vertex weights in Fig. 1 . On the random lattices the corresponding terms in the 
action are :(a) ($t)3 ) ( b ) $3 ^ ($t)2 $) ( d ) $2 $t 
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